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Abstract 

New nonlinear equations for field correlators in the large Nc glu- 
odynamics are derived. Nonperturbative confining correlators are 
shown to satisfy a system of nonlinear scale-invariant equations. Keep- 
ing only the simplest, bilocal correlators one can find the exponential 
large distance behaviour for both independent correlators D and Di 
and express coefficients in the exponentials - the gluonic correlation 
length Tg - through the string tension. Numerically these results are 
in agreement with lattice data. 

1 Introduction 

The main unsolved problem of QCD and of the modern field theory in general 
is the quantitative treatment of nonperturbative phenomena. There is no for- 
malism yet (e.g. a system of equations) in which nonperturbative quantities, 
like masses and condensates in QCD, are computed from the first principles. 
By now there are two main approaches to this problem: i) lattice calculations; 
ii) modelling the nonperturbative vacuum with e.g. quasiclassical solutions. 
In the first approach one can indeed calculate numerically hadron masses 
and even more fundamental quantities like field correlators [1]. However, 
here the physical mechanism creating nonperturbative phenomena is kept 
hidden behind the formidable numerics. The second approach is necessarily 
model dependent and approximate. Also there is no successful model of the 
QCD vacuum, which describes both basic nonperturbative phenomena: con- 
finement and chiral symmetry breaking (CSB) (however, the dyonic model 
recently suggested for QCD [2] might be promising in this respect). 
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The main characteristics of nonperturbative (NP) phenomena in QCD is 
that they create there own mass scale. By this fact the scale invariance is vi- 
olated and through the scale anomaly the NP shift downwards in the vacuum 
energy density is produced [3]. Therefore one may expect that NP solutions 
appears as scale violating solutions of some nonlinear scale-invariant equa- 
tions, while pcrturbativc solutions of the same equations are obtained by 
expanding in powers of coupling constant. 

Recently equations of this type have been written for gauge-invariant 
quark propagators [4], where the kernel of the equations contains the fun- 
damental string. The equations [4] do not violate symmetry (in this case 
the chiral symmetry), while NP solutions break chiral symmetry. In [4] the 
gluon-field correlators were taken as an input and the problem how to find 
the correlators from the first principles was not solved there. Here we are 
proposing selfcoupled equations for field correlators in the case of gluodynam- 
ics. These equations are also of Dyson-Schwinger type but will be written 
for the gauge-invariant amplitudes, which automatically display at large dis- 
tances the formation of adjoint string in the large Nc limit. 

To write these equations we divide all — 1 gluonic fields into two 
groups: a small group of r fields 6^, i ~ 1, ■■■r, which later be called "valence 
gluons", and a large group of fields A^, a — r + 1, ...N^ — 1, which form a 
background for valence gluons. 

As will be shown below this background creates in the confining phase 
after vacuum averaging an adjoint string. Hence, physically the problem 
reduces to the motion of valence gluons in the field of the string. Since the 
string is a white object, the action of the string on valence gluons generates 
color-diagonal mass operator and the resulting equation for the valence gluon 
propagator simplifies considerably. We shall find out that these equations are 
nonlinear and selfcoupled, so that one can look for selfconsistent solutions for 
the valence gluon propagator and the field correlator. Inserting the resulting 
field correlator in the kernel, one realizes the selfcoupling mentioned above. 

Some properties of the solutions can be found even before the actual 
solution of equations e.g. the exponential large distance behaviour of field 
correlators and the coefficient in the exponent - the so-called gluon corre- 
lation length Tg - can be computed through the string tension or gluonic 
condensate. Our number for Tg is in agreement with lattice calculations [1]. 

Thus one finds a solution depending on a single scale. This scale should 
be given as an input to define the NP part of the theory completely. 
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We shall not touch in this letter the perturbative expansion and renor- 
malization of equations obtained here, referring the reader to a future pub- 
hcation. The paper is organized as follows. In section 2 the separation of 
gluonic fields into "valence" and background parts is done and the effective 
action for the valence gluons is written explicitly. In section 3 the mass op- 
erator for valence gluons is found after averaging over background fields. In 
section 4 the connection of field correlators with the valence gluon propagator 
is established and it is used to find the asymptotic behaviour of field corre- 
lators and to compute the gluonic correlation length Tg through the string 
tension or gluonic condensate. In Conclusion a summary of results is given 
and possible extensions are discussed. 

2 Effective action for valence gluons 

We start dividing the gluonic fields A^(a = 1, ...iV^ — 1) into a group of valence 
gluon fields with the color index a running a small set of colours: a — 
1, ...r, while the rest of fields, for which we keep notation ^4°, belong to the 
large set, a = r + 1, ...A*"^ — 1. Number of colors A*"^ — 1 is assumed to be large, 
and r -C N^ — l. Note that this division does not separate perturbative gluons 
from a NP background, as is usually done in the background perturbation 
method [5]. Both 6^ and may contain perturbative and NP components. 
Hence, = A" -|- 6" and one can write for the Euclidean action density 



SE{A) = SE{A + b) = 

= + S^'-\h, A) + S^^'^\h, A) + 5(2'°) (6)+ 

(6, A) + 5(3'^) (6, A) + {b, A) + ^(3) (^,) + ^(4) ^ 



(1) 



where we have defined 



S^^\b,A) = 

5(2'°)(6) = 6;( 



ldl6,, + ld,d.)bl 



(4) 



(2) 



(3) 



S^'''\b,A) = ^gr'% 



[2AldxS^,+ 5^ Al - Ald,]b\ 



k 



V 



(5) 
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S^^\b,A)=gr''%lbtF;,{A) (6) 

S^'''\b, A) = ^x"''''\blblAlA\ - b^blAlAl) (7) 

S^^^'\b,A) = blblx''^HlAl (8) 

S^'\h) = Iv.blx'^^'^^b'X. S^'^ = l9f%b'xdxb'; (9) 
In (2)- (9) the following notations are used 

x'"^''' = 9'Y. nr'"^ (10) 

c=l 

Note also that we have denoted color indices of 6^ fields by i, k,l, ... running 
from 1 to r, and indices a,b,c in Aj^ running from r + 1 to — 1. 

As a next step consider averaging over fields in some expression, like 
the following 

D,.M^, y) =< F;,(a;)$"'^(a;, y)F^Ay) >b,A, l<a,P<N',-l, (11) 
where $ is the parallel transporter in the adjoint representation 

^''"{x, y) = {Pexpig f A^dz^p. K = K^"^ = -^/"'' (12) 
Jy 

In computing the matrix element (11) one can integrate first over fields 
and then as a last step integrate over fields 6^. It is essential that $ contains 
both fields A^j^ and 6^ and we shall choose the so-called modified coordinate 

gauge [6] to express in terms of field strength F^^y. Let us now specify 
this gauge and to this end choose x, y in (11) on the 4-th axis, so that 
X = (0, 0), y = (T, 0). In the chosen gauge we have 

AAz,zA— I daZiF^Aaz, zA, (13) 
Jo 

Ak{z, za) = / adaZiFik{az, Z4) (14) 
Jo 

From (13) it is clear that $(0, T) = 1 and parallel transporters can be omit- 
ted. 
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The basic quantity to be considered below is the gluonic propagator 

G%{x,y)^<l/^ix)b':,iy)>,,A, (15) 

where the averaging over fields b, A should be done with the weight exp{— J SEd^x) . 
The averaging over fields 6, will be done in two steps. First, the averaging 
over An for some operator N can be written as 



< 



N{A) >A= ^ j DA^expi- j d''xS^'^\A))N{A) (16) 



and therefore the averaging over A^, 6^ in (15) implies the definition of the 
effective valence gluon action Seff{b), 



=//W=<e-Ej5'd*->^ (17) 



with i = (1), (20), (21), (F), (2, 2), (3, 1), (3), (4). As a next step the averaging 
over fields b can be writen for some operator T{b) as 

< T{b) i I exp{-S,ff{b))T{b)Dbn (18) 

In calculating Seff{b) in (17) we shall use the cluster expansion theorem and 
neglect in the first approximation all higher cumulants beyond quadratic ones 
(discussion of this approximation will be given below in this section) . Then 
one obtains 

where double angular brackets denote the cumulant: 

< ^2 S^>-<S>'' (20) 

It is easy to show that within our approximation the terms S'^'^'^\ S^^\ S^^^ 
contribute only to the linear average < >^, while S'^^'^\ S^^\ S'^^'^^ - 
only to the quadratic one. 
Hence one has 

S,ff{b) = J (5'(2'0)(6) + 5(3)(6) + 5W(6))d4^ + y < ^(2,2) >^ _ J2 L^"^: 

(21) 
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where L^'''\i = 1; (2, 1); (3, 1); F) are quadratic averages of the corresponding 
terms 5*^*^ and also the term, which will be important in what follows, is 
present: 

-V^{x)ht{x)<Al{x)Al{x)>A} (22) 
Now using (14), (15) one can write 

/ < Al{x)Ai{y) >A= ^Jaa'(^,Z/), (23) 

where Cl — ^^n'^ important function J^,^(a;, y) is introduced: 

Jxx'{x, y) = C Xidaa^^^ C dPy^P^^'^g^^-^ < Fa{xa)Fkx'{y(3) >a (24) 

Jo Jo I\c 

Here a^^\ /3^^^ = 1 for A = 4 and = a,/?^ = /3 for A = 1,2,3. Below 
we shall use D and Di correlation functions defined through the following 
decomposition [7]: 

^ < F,x{u)Fkxiu') >A^ {SikSxx' - SiX'Skx)D{u, u') + Ag^^iu, u') (25) 

with the factor A*^^) proportional to the full derivative and therefore not 
contributing to the confinement and to the string creation; we shall disregard 
it below. 

Now we can also rewrite (22) for < 5"^^'^^ > through the kernel 7^,^: 

< S^'''^ >A= -^^hl{x%{x){5,,Jxx{x,x) - J^Ax,x)). (26) 

At this point it is necessary to argue why in Seffip) (21) we have omitted 
the higher-order cumulants and shall also omit the term L*^^'^-*, quadratic in 
5'(2,2)_ fjj^g reason is that those terms contain quartic (and higher) terms of 
the type ^ ^4^ » which due to (13), (14) are expressed through <^ 
As it was argued in [8] these quartic and higher cumulants do not contribute 
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significantly to the string tension which can be written (for the fundamental 
string) as 



For static charges of higher SU(3) representations the correlator D{x) is 
proportional to the quadratic Casimir operator C2 [8] in good agreement 
with lattice data [9]. The contribution 0(-C F"^ ':$>), which produces the (6*2)^ 
dependence in the string tension, has not been seen in the static potential 
of adjoint charges on the lattice [9]. Additional arguments come from the 
string profile calculations [10] where using only D{x) one reproduces lattice 
data with good accuracy. Hence, we expect that the contribution of higher 
cumulants is suppressed at least in those two instances, whereas confinement 
(nonzero a) is already present due to the Gaussian correlator D{x). 

Therefore to obtain the selfcoupled equations for the lowest cumulants 
we confine ourselves to functions D{x) and A^^^ while at the next stage 
one can express triple and quartic correlators through the lowest ones and 
so on. In this way one gets an expansion in powers of the parameter p = 

{g\J< Fff^ >Tg), where indices i, k are fixed and refer to the surface of inte- 
gration, as in (27). Inserting phenomenological values of gluonic condensate 
[3] and the lattice value of Tg [1], one obtains p ~ 0.1 which implies a fast 
convergence of cluster expansion. Note that the case of instantons, where 
the cumulant series is not converging, requires a special treatment, which 
was done e.g. in [11]. 

3 Mass operator and nonlinear equation 

The effective action Seff{b) in equation (21) contains terms which are quadratic 
in b{i.e.L^^\ < S'*-^'^^ >, S**-^'*^^), and in addition also cubic and quartic terms 
in 5''^'^^ 5''-^-* and L^'^'^\ L^^\ and the sixth power term L'^^'^\ Keeping aside 
5'(3) c;(4) other terms and write the effective mass term 

M in the same way as it is usually done in the Dyson-Schwinger formalism 
[12]. Writing diagrams based on the vertices L(^'^\ L^^^ and L^^'^^ one ob- 
tains (in the leading order of large N^) the Hartree-Fock approximation for 
M generated by those vertices, which are obtained doing a Wick pairing for 
every extra pair of operators with replacement b^b^, — > G^i, (15). 




(27) 
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In the same approximation the valence gluon Green's function Gf^u (16) 
is diagonal in color (since the resulting mass operator M will be shown to be 
is diagonal) and one can rewrite (15) as 

< 6^(x)6^(y) >,^A= SapG,,{x,y). (28) 

As a result one can represent M in the following form, 

M - M(2'2) + m(2'^) + M(^) + M(3'^) + M^^^ + AM, (29) 

where M^'^ correspond to 5'^*^ and e.g. 

Mj^fHx,y) = ^S^^){x-y)[Jxxix,y)d,, - J,.{x,y)], (30) 
Mj^/\x, y) = -Ajxx'ix, y)Gp,ix, y)Nxp^{x)Nx'aM (31) 

with 

T 

and the parameter ^ = — — , (32) 

N„ — 1 



c 



and 



M'^Kx.y) - < D'^FUx)Df,{y)Fi,M >A- 



c 

1 a (9 



= -iV|3I&l9^ < >. +0« FFF » (33) 

Here in (32) we have omitted higher cumulant contribution and all other 
M^*) proportional to ^ will be disregarded in the large Nc limit. 

The terms S^^\ S^^^ generate their own contribution to M^^,, which we 
denote AM. We stress now that all terms in M containing the valence 
Green's function G^i,, acquire a factor ^ for each G, the same also refers 
to AM. Except for ^ the dependence in all terms in M cancels and all 
terms are 0(iV°). Therefore in the first approximation one can neglect all 
terms containing ^ and is left with 

M = M(2'2) + MW + 0(l/iV2) (34) 
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Then Dyson-Schwingcr equation for the valence gluon Green's function can 
be derived from (21 ) in a usual way [12]: 

{-dXp + d^dp)Gp,ix,y)+ J M^p{x,z)Gp,{z,y)d'z ^ 5^^\x -y), (35) 

here M is given in (24-33) or at large Nc approximately in (34). 

Equation (35) is the central point of the paper. In the rest part of it 
we analyze equation (35) in more detail and express field correlators D, Di 
through G/^i,, thus closing the set of equations. 



4 Analysis of the selfcoupled equation (43) 
and simple estimates 

To compute D and A^^\ as defined in (25), through G^i, it is convenient to 
choose, in each term of the sum over color indices in (11), the color indices 
a, (3 as belonging to the small set of valence colors. Hence fields F^i^ will 
consist of bn only and (11) can be rewritten as 

d 

+perm + gf'' < - — b'^{x)^'^''{x,y)bl{y)bl{y) > +perm+ 
ox^ 



+ g'r'^r' < bi{x)bt{x)bl{y)b'M >} = ^ \^ ^ {h + h + h] (36) 



By perm in (36) we denote the sum of terms which obtain from the preceding 
one using the rule 

(liuXa) + perm = +{iJ,u\a) — {uixXa) — {ixuaX) + {ujiaX) (37) 

In (36) it is also assumed that finally one can put x = (00), t/ = (2^,0) and 
therefore ^{x,y) in the gauge (13) can be replaced by unity. 
The r.h.s. of (36) can be written symbolically as 

di,d^G^a(x,y) + G^^{x,y)G^a{x,y) +perm. + 0(< FFF >) + 0(/). 
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Since the asymptotics of G{x, y) at large |a;4 — 1/4I ( as is seen from the integral 
equation (35) with the kernel (34)) is exponential, the same is true for the 
l.h.s. of (36), i.e. for D(x) and A^-^^ This fact is in agreement with lattice 
data [1]. 

The equation (35) and the kernel M do not contain dimensional param- 
eters and arc scale-invariant. Therefore going from large distances, where 
nonperturbative contributions are dominant, the scale appears spontaneously 
in the nonperturbative solution. All other nonperturbative paramerters, like 
gluonic condensate ~ (-D(O) + -Di(O)) are expressed through this prescribed 
scale parameter. 

If one, however, starts from small distances and develops the perturbative 
series like 

G^,{x, y) = G^o) - G(°)MG(°) + G(°)MG(°)MG(°) - ... (38) 

where G^°^ is the free solution of (35), then the perturbative divergencies 
require regularization and the renormalization, which introduce their own 
perturbative scale - normalization mass ji or Aqqd ■ Fixing this scale, one 
can extrapolate to larger distances and find nonperturbative solution defined 
at this scale. We shall discuss this procedure in detail in a subsequent publi- 
cation, and now look at the properties of the pure nonperturbative solutions. 

We shall be interested in the large time behaviour of gluon Green's func- 
tion and to make an estimate we first keep in the kernel M (34) only the local 
part since it grows at large distances, while M^^^ tends to a constant. 

We assume at this point (and later confirm it in calculations) that D{u) 
has an exponential form 

D{u) = Doexp{-u}\u\), u = l/Tg. (39) 

Prom (27) one has a — and J\\ is 

J,,(x,x) = ^((1)s + (|)m) (40) 

We have used in (40) the subscripts -E, M to specify the contribution from 
colorelectric and colormagnctic fields respectively. The latter is gauge- de- 
pendent if one keeps only Gaussian correlators and neglects all higher-order 
ones. 
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One can diminish this gauge dependence taking into account that the 
inclusion of all correlators creates for every Wilson loop the surface of the 
minimal area. In our case of the large T asymptotics this surface, built 
on the trajectory of the valence gluon and the parallel transporter along the 
time axis, is time-like. This means that (unless large angular momenta J are 
considered) the spacial projections Yl,ij of the minimal surface arc not growing 
with time T, and therefore magnetic contributions to the linear confinement 
term in (40) should be cancelled in the sum of all correlators. Hence, the 
choosing of the minimal area surface in the path-integral formalism [13] is 
equivalent to the temporal surface gauge [14] in the formalism presented here 
and in [4]. As a practical outcome the magnetic term (|)m in (40) vanishes 
in the temporal surface gauge. 

Prom the structure of (35) one can deduce that G^y{x, y) can be expanded 
in a complete set of transverse gluon states 

G^.(^,Z/)=EC«n(xX(Z/) (41) 
n 

where eigenvectors Un{x) satisfy an equation 

- dlun{x) + ^Jxx{x, x)un{cx) = 0. (42) 

Writing Un{x) as 

Un{x) = exp{-LLlnX4)(Pn{x) (43) 

one obtains equation for (fn{x) 

i^l^n = (-^ + C\x\)^n (44) 



with 



dx^ 



C ^ (45) 



'2 

The analysis of the equation (44) was done in [15] with the result 

ujI = a^{Cf'\ (46) 

where a„ are numerically computed numbers, for the lowest state n = one 
has ttQ = 2.33... 
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At this point it is important to use the consequence of selfcouphng in (35) 
and (36). Namely, as follows from (41) and (43), the propagator G^i, decays 
exponentially at large time T. As a consequence of (36) and the analysis 
done above, the function D{T) decays at large T with the same exponent, 
i.e. 

uj = ujq (47) 
Inserting (47) into (46) one obtains 

u = {222>fl\l^ (48) 

V ZTT 

and with a = 0.2GeV^, the final result is 

u^l/Tg^ IMGeV (49) 

This value agrees well with the lattice measured value Tg = 0.22/m(a; = 
OmCeV) [1]. 

Thus we have demonstrated in this example that different parameters of 
the nonperturbative correlators can be expressed through each other via a 
dynamical computation. 



5 Conclusions 

The analysis of the previous sections was intended to show that selfcoupled 
equations (30), (34-36) indeed contain a nontrivial dynamics, which allows 
to obtain nonperturbative gluon propagators and correlators. Numerical so- 
lution of these equations is possible and is now in progress, in this paper 
however we confine ourselves to simple estimates and qualitative analysis. 
The characteristic behaviour of the gluon propagator G{xi,y4,x,y), which 
one finds from (35) and (43), is exponential decay in — y^l and the faster 
than exponential decay at large |a;| or \y\. The latter property signals the 
confinement of the gluon, since the adjoint string at large Nc does not allow 
the gluon to go far from the origin. The parallel transporter along the tem- 
poral axis, present in our definition of the gauge-invariant gluon propagator, 
is actually a signal of the static adjoint charge at the origin, therefore the 
gauge invariant quantity Gfj_i, is in fact the propagator of the white system 
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made of the gluon and the static charge. Hence the mass uj in the exponential 
decay of G^J,^, in time is equal to the mass of this white system. 

Connection between G^^^ix/^^ y^, 00) and D{xi — y^ and I?i(x4 — ^4) found 
in section 4, tells us immediately that D and Di decay with the same domi- 
nant exponent with mass uo. 

In the subsequent paper we shall discuss the perturbative contents of our 
equations (35) and of the effective action S^ffih) (21). 

One of the most important problems to be solved in the framework of 
the suggested formalism is the problem of the deconfining phase transition. 
The latter may occur or due to the increase of temperature or due to the 
addition of scalars with the Higgs-type interaction. In the last case one 
can get a system of equations of the type of (35) for the gluon and scalar 
propagators which describe two possible phases: confinement and deconfined 
massive gluons. We plan to analyze these questions in forthcoming papers. 

The author is grateful to A.Di Giacomo, H.G.Dosch, M.Schmidt and 
V.l.Shevchenko for useful discussions and to A.M.Badalian for valuable re- 
marks. 
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